arXiv:1210.2027vl [physics.optics] 7 Oct 2012 


Effective spontaneous VT -symmetry breaking in hybridized metamaterials 


Ming Kang 1 , Fu Liu 1 , Jensen Li l 2 0 

1 Department of Physics and Materials Science, City University of Hong Kong, Tat Chee Avenue, Kowloon Tong, Hong Kong, China 
2 State Key Laboratory of Millimeter Waves, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong, China 
(Dated: October 9, 2012) 

We show that metamaterials can be used as a testing ground to investigate spontaneous symmetry breaking 
associated with non-Hermitian quantum systems. By exploring the interplay between near-field dipolar coupling 
and material absorption or gain, we demonstrate various spontaneous breaking processes of the PT -symmetry 
for a series of effective Hamiltonians associated to the scattering matrix. By tuning the coupling parameter, 
coherent perfect absorption, laser action and gain-induced complete reflection (n reflector) by using an ultra- 
thin metamaterial can be obtained. Moreover, an ideal PT -symmetry can be established effectively in a passive 
system by using metamaterials. 

PACS numbers: 78.67.Pt,42.25.Bs,78.20.-e 


Metamaterials, periodic arrangements of artificial subwave¬ 
length atoms, have paved the way for realizing different novel 
concepts and applications in optics, such as negative re¬ 
fraction, sub-diffraction imaging and invisibility cloakinglQ]- 
0 ], Instead of using the macroscopic material parameters 
for description, many of the optical properties, including 
metamaterial-based chirality la lal and electromagnetically in¬ 
duced transparency (EIT)i7l4l411. can be understood by study¬ 
ing the coupling between different magnetic and electric artifi¬ 
cial atoms in the so called hybridization picture J3,[l2lE,0]- 
Similar to optical lattices with controllable counhng |[l7l - [l9ll . 
the design flexibility in the coupling between artificial atoms 
opens up a new perspective to realize quantum simulators 
to achieve effects which are difficult to obtain in quantum 
svstems lfl7l[l8ll . Here, we theoretically investigate the pos¬ 
sibility in constructing a metamaterial system with an effec¬ 
tive PT -symmetric Hamiltonian which can still show real 
eigenenergies. One of the most interesting properties of a 
PT -symmetric Hamilton is the existence of a phase tran¬ 
sition arising from spontaneous symmetry breaking by tun¬ 
ing an external parameter, defined as the coalescence of two 
eigenvalues and their corresponding eigenvectors lEfill . Un¬ 
like the previously considered symmetry breaking effects of 
metamaterials ll2lll . such a spontaneous symmetry breaking is 
implicit in nature. It means that the phase transition occurs 
without structurally destroying a specific symmetry (e.g. mir¬ 
ror, rotational symmetry). The external parameter being tuned 
still respect the PT -symmetry. 

During the past few years, f’T'-symmetric Hamiltoni¬ 
ans have in fact been extensively investigated within the 
framework of optics. The most direct analog is repre¬ 
sented by a pair of optical waveguides with a balanced gain 
and loss profile ll9t 122143 1 II . Spontaneous PT -symmetry 
breaking has also been theoretically and experimentally 
investigatedjjjl Izj-fiTtl. Intriguing phenomena, such as 
asymmetric transmissi onlll9[El . power oscillation |0], uni¬ 
directional invisibility 1191 |29[], coherent perfect absorption 
and laser lf3oh have been predicted. Moreover, extending from 
an ideal PT -symmetric system with balanced gain and loss, 
non-ideal PT -phase transition of a passive system has also 


been proposed and experimentally demonstrated E 0. 
Unconventional loss-induced optical transparency is demon¬ 
strated around the phase transition point ll25n . Actually, PT- 
symmetry has not been discussed within the context of meta¬ 
materials. By employing metamaterials, the larger flexibility 
in designing the coupling between different artificial atoms 
gives us additional possibilities in realizing non-trivial PT- 
symmetry. Here, we show that by mapping a metamaterial 
consisting of a bright electric and a dark magnetic atom to 
an effective Hamiltonian, we can establish an ideal sponta¬ 
neous symmetry breaking transiting from little absorption to 
complete absorption even in a passive system. It relies on the 
inverse relationship that scattering is much larger than loss for 
a bright electric atom while the reverse is true for a dark atom. 
Such a system has been previously considered to achieve EIT. 
However, we are interested in the phase transition to obtain 
complete absorption llE lEl across a threshold in the cou¬ 
pling parameter, in contrast to the high transmission for EIT 
emerging from the strong coupling regime iQ]- Moreover, by 
incorporating a small material gain into the dark atom of the 
ultrathin metamaterial, we can obtain laser action lE - tilfill and 
gain-assisted complete reflection by considering the sponta¬ 
neous symmetry breaking effect of a class of Hamiltonians 
associated with the scattering matrix. 

The discussion of PT -symmetric Hamiltonians originates 
from Bender and colleague’s work, showing that a Hamilto¬ 
nian H respecting a parity-time symmetry, i.e. [PT,\1\ = 0, 
can have completely real eigenenergies ll20ll . To find the ana¬ 
log of PT -symmetry in a metamaterial system, the simplest 
way is to look for a pair of coupled artificial atoms, since 
their response can be understood within a hybridization pic¬ 
ture in quantum mechanics as a two state problem, and the 
PT -symmetry in two-state quantum system is already well 
studied lEflll . Here, we consider a metamaterial constituting a 
planar array of atoms. Each unit cell consists of an electric 
atom coupled to a magnetic atom through the near fields. For 
simplicity, we assume normal incidence and there is no polar¬ 
ization conversion for the incoming plane waves. The optical 
response of such a metamaterial can be expressed using the 
scattering matrix S as 


( 1 ) 


where a± and b± are the E-field amplitudes of the incoming 
and outgoing waves in the positive and negative directions. 
The same scattering process has a microscopic representation 
using the atomic responses and their associated far-field ra¬ 
diations. By further assuming each atom is dominated by 
only one dipolar resonating mode, the atomic responses can 
be summarized as 
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while the outgoing fields and the atomic responses are related 
to each other by 
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In the above equations, y, = yf + y l ° ss (i = 1,2) represents 
the total resonating linewidth for each atom where yj, a pos¬ 
itive number, represents the part of scattering and y , ° ss is the 
part due to absorption loss when it is positive. We also nor¬ 
malize the electric dipole p and the magnetic dipole m by the 
oscillator strengths as in Ref.7, and assume the two resonant 
frequencies at the same frequency (Oi, res = <oq. In fact, Eq. (4) 
defines a series of Hamiltonian H„ with meaning in the fol¬ 
lowing. By taking an appropriate arithmetic average of Eq.(2) 
and Eq.(3), one can prove that det(H n - col) = 0 yields an 
eigenvalue (n - 2) In of matrix S. We can also obtain 


det(H 2 - oA) _ ^ det(Hi - oA) 
det(Ho - oA) det(Ho - oA) 


(5) 


For example, det(Ho - a A) = 0 leads to the pole of S, which 
corresponds to laser action |[30l l34l - [36ll . det(H 2 - oA) = 0 
implies a zero of S to completely absorb the incoming EM 
waves, namely coherent perfect absorption fCPA) l3oi [ 32 L[33ll . 
For a metamaterial with mirror symmetry in the propagation 
direction, the two eigenvectors of S are (a+, «_) = (1,1) with 
eigenvalue t+r and (1,-1) with eigenvalue t-r by defining t/r 
as the complex transmission/reflection coefficient for electric 
field with respect to the mirror plane. If we further confine our 
discussion to the coupling between a bright electric and a dark 
magnetic atom, we have the additional requirement t = 1 + r. 
Then, det(Hi - oA) = 0 implies a complete reflection r = -1 
and t = 0. In general, a particular ^/eigenvalue of S denotes 
a particular optical phenomenon and it corresponds to solv¬ 
ing the eigenvalue problem of H n . If the eigenvalue of H n 
is real, the phenomenon can occur at a real frequency with a 
pole/zero on the corresponding optical spectrum. Therefore, 
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by considering the spontaneous symmetry breaking process 
of a VT -symmetric H n , one can observe the phase transi¬ 
tion when the real eigenfrequencies are turned into complex 
conjugates or vice versa. For our 2x2 Hamiltonian, the P- 
operation switches rows and columns of the matrix while the 
T-operation takes complex conjugate on the elements. A VT - 
symmetric H„ is defined by [VT, H n ] = 0. We have to empha¬ 
size that VT -symmetry being discussed here is an effective 
one through the establishment of an effective VT -symmetric 
Hamiltonian, which is different to the approach employing 
complex notentials lEoll . 

As an example, a VT - symmetric H 2 corresponds to the 
condition y 2 - 2y 2 = 2yj - yj, which can be satisfied by a 
representative hybridized metamaterial with a bright electric 
atom (scatters most) and a dark magnetic atom(absorbs most). 
For a fixed set of parameters, e.g. j\ = 15 THz, y l ° ss = 1.75 
THz (yj = 13.25 THz), y 2 » y 2 ss , this condition is repre¬ 
sented by the vertical red line in Fig. 1 at y 2 = 11.5 THz 
to indicate a VT -symmetric H 2 . Then, the spontaneous sym¬ 
metry breaking process can be observed by changing the cou¬ 
pling coefficient k (still respecting the VT -symmetry) across a 
threshold value. The eigenfrequencies of VT -symmetric H 2 
are real(solid red line) when k > y 2 , and becomes a com¬ 
plex conjugate pair(dotted red line) when k < y 2 . The split¬ 
ting of the real/imaginary part of the eigenfrequencies with 
varying k is showed in the inset. This behavior reveals itself 
through the absorption coefficient when the metamaterial is 
excited by the incident waves of (a+, a_) = (1,1), giving CPA 
for k > y 2 at the two real frequencies. We note that such 
a configuration of coupling between a bright and a dark atom 
(without the VT -symmetric condition satisfied in general) has 
been used to give EIT in the strong coupling regimefl. We 
endow the same system to visualize effective VT -symmetry 
breaking and focus on its optical properties in the regime near 
the critical point. In the following, we demonstrate the inter¬ 
esting optical properties associated with VT -symmetry break¬ 
ing with metamaterial structures. In particular, undergoing the 
phase transition, the metamaterial becomes perfectly absorp¬ 
tive with real frequencies. The bright atom (inset of Fig. 2(a)) 
is a silver strip of length LI = 150 nm, width W 1 = 50 nm 
and thickness t = 20 nm. Fig. 2 (a) shows the electric field 
near one end of the strip with numerical simulation (symbols) 
and theoretical model (lines). The numerical simulation is car¬ 
ried out using finite element analysis (Comsol MultiPhysics) 
with Drude model of silver being e = &o - co 2 p {(o 2 + iy p io)~ l , 
where co p = 2.1961 X 10 15 Hz, y p = 4.3439 X 10 12 Hz, and 
6oo = 3.7. The theoretical result is fitted using a Lorentz model 
E s oc (cjo -co- iyi) _ 1 (y 1 - Y° ss )Eq, where u>o = 402.5 THz, 
yi = 15 THz, and y( oss = 1.75 THz, E v (0) is the scatter- 
ing(incident) E field. y° ss can be estimated by the probing the 
response of the bright atom with and without metal loss. The 
dark atom is a magnetic resonator tuned to the same resonat¬ 
ing frequency with two parallel silver strips and a dielectric 
gap material (e = 2.25 + ie" with e" introduced to represent a 
slight absorption or gain) in between. The geometric param¬ 
eters are L2 = 100 nm, W2 = 30 nm, t = 20 nm, as showed 
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FIG. 1: (Color online) Parametric evolution of the metamaterial for 
different H„(;z = 0,1,2); Inset: Splitting of eigenvalues for H 2 with 
varying k. Vertical solid lines indicate the real eigenvalues case for 
ideal PT -symmetric H n above the phase transition threshold, short 
dotted lines illustrate the complex conjugate pair eigenvalues case for 
ideal PT -symmetric H„ below the phase transition threshold. Inset 
gives the eigenvalues with varying k for PT -symmetric H 2 . Curved 
solid lines show the real eigenvalue condition indicating the optical 
properties for non-ideal PT -symmetric H„, short dashed lines show 
branching condition for non-ideal PT -symmetric H n . 


in the inset of Fig. 2(b). Due to the large quality factor of 
the dark mode resonance, a small e" of the gap material can 
effectively alter the half-linewidth y 2 . As shown in Fig. 2 (b), 
numerically calculated field (symbols) near one of the metal 
strip ends is again fitted well with theoretical Lorentz model 
(solid line) for e" = 0,-0.1 (hollow, solid symbols). Com¬ 
paring between numerical and analytical results with different 
e", a simple approximated formula between y 2 and e" can be 
established as y 2 = 1.75 + 22,5e". Now, we hybridize these 
two atoms within a metamaterial with a square lattice of pe¬ 
riodicity d = 400 nm and a distance h in separating the two 
atoms, as shown in Fig. 2(c). The change of distance h be¬ 
tween the two meta-atoms can tune the coupling k between 
them, but there exists no simple relationship between them, as 
shown in Fig. 2(d). The k can be calculated from the bright 
meta-atom probe response for e" = 0 by changing the dis¬ 
tance h between them. The fitted formula can represent main 
variation character of k, but exhibits slight deviation in large k 
region. 

Now, we illustrate the spontaneous symmetry breaking with 
fixed geometry of the meta-atoms with constant a>o, yi, and 
y/ojj The variable parameters are the coupling coefficient k 
(or equivalently h ) between the two atoms and y 2 (through 
the tuning of e"). We can then choose a particular y 2 to sat¬ 
isfy PT -symmetry requirement of H n . For a PT -symmetric 
H 2 , we need y 2 = 2 y[ - yi = 11.5 THz, i.e. e” = 0.433, 
which is a passive configuration to display the effective PT 
symmetry. In the fight of the above theoretical discussion, the 


FIG. 2: (Color online) Analytical (solid line) and numerical (scatter 
line) spectral response of the E-field probe for bright(a) and dark(b) 
meta-atom to establish the atomic Lorentz model with fitted param¬ 
eters (cu 0 THz, yi THz, yand y 2 ). (c) and (d): The hybridized 
metamaterial under consideration and the relationship between the 
coupling coefficient k and the distance h between the two meta- 



FIG. 3: (Color online) Observation of effective spontaneous PT- 
symmetry breaking for H 2 and H| Analytical model (with fitted 
parameters from Fig. 2) (a) and numerical(b) outgoing \t + r\ 2 spectra 
versus frequency and coupling distance. Analytical(c) and numeri- 
cal(d) transmission spectra versus frequency and coupling distance. 


eigenvalue problem for a PT -symmetric H 2 is equivalent to 
the zeros of the S, which can be easily observed in the co¬ 
herent outgoing \r + t\ 2 spectra when we change the param¬ 
eter h. The PT -phase transition point in this case occurs at 
Kpj = 11.5 THz (h = 46.0 nm). Below the critical point, 
kpt < 11.5 THz (h > 46.0 nm), the eigenvalues of H 2 is 
a complex conjugate pair, there are no crossing in real fre¬ 
quency axis, so we cannot observe CPA in this region, shown 
in Fig. 3(a). As we increase (decrease) k(H), at kpt = 11.5 
THz(/z = 46.0 nm), the two eigenvalues of H 2 collapse to 
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one point in the real frequency axis, i.e. CPA(|r + r| = 0) 
emerges in the spectra. When k passes over the kpt, there 
are two real frequencies satisfies the eigenvalue solution of 
H 2 , CPA branching phenomenon can be clearly shown in Fig. 
3(a) from the analytical results obtained from Eq. (2) and (3) 
with the atomic model established from the previous simula¬ 
tions shown in Fig. 2. The same phenomenon can be ob¬ 
served from the numerical simulations (with microstructures 
included) shown in Fig. 3(b). Real parts (white dashed line in 
Fig. 3(a),(b)) of the eigenvalues of H 2 agree well with the CPA 
positions in the spectra. In both Fig. 3(a) and Fig. 3(b), the 
absorption loss (1 - \t + r \ 2 for symmetric incoming waves) of 
the hybridized system increases dramatically with increasing 
coupling (decreasing separation) below kpt- At kpt, the pres¬ 
ence of m interferes with p to the right size and gives r = 
and t = \. This is the onset of CPA phase. Going beyond 
the transition point (larger coupling), the CPA peak splits into 
two. On the other hand, this splitting makes both atoms can¬ 
not be efficiently excited at to o and loss decreases at that fre¬ 
quency. In fact, the non-monotonic function of loss against 
coupling coefficient is a direct manifestation of spontaneous 
symmetry breaking, which is similar to the loss-induced opti¬ 
cal transparency in a passive FT -symmetric system deviated 
from the ideal FT -symmetry J3- Here, using metamaterials, 
we can actually establish an ideal FT -symmetry in a com¬ 
pletely passive system. CPA becomes possible besides the 
counterintuitive trend of loss function. 

Besides H 2 in revealing FT -symmetry breaking process, 
Hi is also a potential candidate to display the same process. 
As predicted in Eq. (5) and Fig. 1, by taking into account 
the eigenvalue problem of Hi, we need to fulfill the case 
j 2 = -y,"” = -1.75 THz (the black vertical line in Fig. 1), 
i.e. e" = -0.156, where potential amplifying media can be 
available in nature ll34l - [36ll . The emergence of effective FT- 
symmetry breaking process can be observed though transmis¬ 
sion spectra, and t = 0(r = -1) in spectra implies real fre¬ 
quency solution for eigenvalue of Hi. It indicates the system 
without absorption loss, A = 1 - \t + r| 2 = 0, and transmis¬ 
sion, t = 0, with perfect reflection with and a w-phase shift, 
r = — 1. The phase transition point in this configuration is at 
kpt = 1-75 THz (h = 98.4 nm). Illustrated in Fig. 3(c), and 
numerically verified in Fig. 3(d), the phase transition point 
is at h = 97.5nm, t = 0 (r = -1) in the transmission spec¬ 
tra agrees well with the real parts of the eigenvalues of Hi. 
Below the phase transition point (k < xpr), the eigenvalue of 
Hi implies no cross point in the real frequency axis, a com¬ 
plex conjugate pair, so transmission is dominated by the bright 
atom dipole oscillation, although t » 0, the eigenvalue of S is 
not -1 (Absorption loss is not zero). Above the phase transi¬ 
tion point k<pt, the eigenvalue of H\ implies two real eigen- 
frequencies, so t = 0 (r = -1) can be observed at two differ¬ 
ent frequencies in the transmission spectra. Real parts (white 
dashed line in Fig. 3(c),(d)) of the eigenvalues of Hi match 
the t = 0 positions well in spectra. The effective spontaneous 
FT -symmetry breaking for Hi can thus be observed in the 
spectra. It is worth to note that the transmission exhibits maxi¬ 


mum value at the position k s = y-y 2 yi = 5.12 THz (h = 73.8 
nm for analytical prediction and h = 72.5 nm for numerical 
results) at ojq, because gain behavior of the system appears, 
when pass over the phase transition point k<pt ■ As we in¬ 
crease k , we can encounter a real eigenvalue of Ho, indicating 
pole of S, which implies the system shows laser action at this 
parameter position lbd 134143611 . as indicted in the cross point 
between the black vertical line indicating FT -symmetric H\ 
and the blue curve displaying real eigenvalue of Ho in Fig. 1. 
The laser action results from non-ideal FT -symmetry of Hq 
whose origin will be discussed in the following. 

For Ho, ideal FT -symmetry condition y 2 = -j\ implies 
that huge gain is required to compensate metal losses. This 
is difficult to be satisfied using natural materials. However, 
in a practical (non-ideal) system, even though ideal condi¬ 
tion is not satisfied, the eigenvalues still split in the same way 
as an ideal FT -symmetric Hamiltonian. Only an additional 
constant bias in the imaginary part of the eigenvalues occurs. 
Similar strategies have been employed in studying non-indeal 
FT -symmetry in passive svstems ll25ll27ll . Then, phase tran¬ 
sition (branching) point is not on the real frequency axis, and 
is shifted into the lower-half of the complex cu-plane. One of 
the eigenvalues (the one with increasing imaginary part) will 
cross the real frequency axis, corresponding to the original 
optical phenomenon, e.g. laser action for the case of Ho as a 
pole of S at k s = y— yi vo Ea l34l - [36h . It is the solid blue curve 
in Fig. 1 and is lying below branching condition (branching 
occurs when k > k<pt) where k<pt = (yi - y 2 )/2(dashed green 
line). In contrary, ideal FT -symmetric Ho has the poles of S 
on the real frequency axis beyond branching condition (ver¬ 
tical solid blue line) at two different frequencies. As an ex¬ 
ample, y 2 = -0.5 THz (e" = -0.1), Fig. 4 gives the outgo¬ 
ing intensity of the system with symmetric incoming light as 
a function of k(K). Faser action collapses to one position at 
oj o, h = 89.4 nm for theoretical prediction and h = 90.0 nm 
for numerical results. The branching phenomenon still occurs 
(without real poles of S) and agrees well with the real parts 
of the eigenvalues of Ho (black dashed lines), as displayed 
in Fig. 4. Faser action and branch phenomenon is the direct 
demonstration of non-ideal P7~-symmetry for Ho. Besides 
Ho, non-ideal FT -symmetry can also occur for Hi and H 2 . 
In Fig. 1, the real eigenvalue of H 2 (zero of S), i.e. CPA, 
is at k 0 = ^y 2 (yl - 2y*° ss ) (solid red curved line) while the 
branching condition is at Kj*r = (yi + 72 - y l ° ss )/ 2 (dashed 
purple line). For Hi, perfect reflection, r = -1, t = 0 (one S 
eigenvalue is -1) is at x_i = ^j— y 2 yj” s (solid black curved 
line) while the branching condition is at K<pr = (y l ° ss - y 2 )/2 
(dashed black line). So, for a non-ideal P7 - -symmetric sys¬ 
tem, the optical phenomenon (e.g. laser, CPA or ^--reflector) 
associated to H„ collapses to a single point (with a k a bit less 
than the one triggers branching) instead of the two branches. 

In the above discussion, k is the external parameter in 
demonstrating FT -symmetry. If we change y 2 with a fixed 
k instead (travelling in a horizontal direction in Fig.l), the 
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FIG. 4: (Color online) Output intensity (color map) of the metama¬ 
terial log 10 \t + r\ 2 revealing spontaneous symmetry breaking of non¬ 
ideal PT -symmetry for H 0 and the real part of the eigenfrequencies 
for H 0 (dashed lines), (a) Analytical model; (b)Numerical simula¬ 
tions with metamaterial microstructures. 



FIG. 5: (Color online) Analytical (solid line) and numerical (dashed 
line) outgoing spectra, as a function of y 2 for the metamaterial at a» 0 . 

metamaterial can function as a laser, the cross point for Ho, 
or a CPA, the cross point for H2 at too. As an example for 
demonstration, h = 75 nm, besides appearance of laser ac¬ 
tion and CPA in spectra, schematically shown in Fig. 5, above 
CPA position, an increase of loss (increasing y 2 ) leads to more 
output of the system, which is similar to loss induced optical 
transparency in a non-ideal passive PT system lEHll . Similar 
counter-intuitive trend also occurs beyond the laser position, 
an increase in gain (decreasing negative y 2 ) causes less output 
observed in the spectra. 

In conclusion, we investigated the effective PT -symmetric 
Hamiltonians using a metamaterial system to reveal sponta¬ 
neous symmetry breaking processes. Theoretical and numer¬ 
ical results clearly confirm the analog of spontaneous break¬ 
ing of PT -symmetry using metamaterials. Our results not 
only open up a new perspective towards testing the intrigu¬ 
ing properties of spontaneous PT -symmetry breaking in non- 
Hermitian quantum systems, but also provide feasibility for 
future applications in plasmonic and metamaterial systems, 
such as laser action ll30U34l - l36l1 . and CPA with sub-wavelength 
thickness ll32ll . Moreover, by using metamaterials, it becomes 
possible to establish an ideal PT -symmetry by using a passive 
system. 
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